function Sal = BeamFormingBruto(Ent,pitch,fs,Zi)



[Nz,Nx]=size(Ent);
Cp = 1.54;
Z = Zi + Cp*[0:1/fs:(Nz-1)/fs]/2;
X = 0:pitch:(Nx-1)*pitch;
Sal = zeros(Nz,Nx);
for i=1:Nz
    for j=1:Nx
        for k=1:Nx
            Dr = sqrt(Z(i)^2+(X(k)-X(j))^2)-Z(i);
            Dt = 2*Dr/Cp;
            Np = round(Dt*fs);
            if (i+Np) <= Nz
                Sal(i,j) = Sal(i,j) + Ent(i+Np,k);
            end
        end
    end
end